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Abstract 

We extract the dynamics imphcit in an algebraic fitted model Hamiltonian for 
the deuterium chromophore's vibrational motion in the molecule CDBrClF. The orig- 
inal model has 4 degrees of freedom, three positions and one representing interbond 
couplings. A conserved polyad allows in a semiclassical approach the reduction to 3 
degrees of freedom. For most quantum states we can identify the underlying motion 
that when quantized gives the said state. Most of the classifications, identifications 
and assignments are done by visual inspection of the already available wave function 
semiclassically transformed from the number representation to a representation on the 
reduced dimension toroidal configuration space corresponding to the classical action 
and angle variables. The concentration of the wave function density to lower dimen- 
sional subsets centered on idealized simple lower dimensional organizing structures and 
the behavior of the phase along such organizing centers already reveals the atomic mo- 
tion. Extremely little computational work is needed. 
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1 Introduction 



In recent years we have developed methods to investigate algebraic models ( spectro- 
scopic Hamiltonians ) for the vibrations of molecules ^ |21 El S El • The Hamiltonians 
reproduce and encode by their construction the experimental data [HI [71 18. .9. lOt fTTl [T^ . 
So far we have treated examples where the existence of conserved quantities allowed 
the reduction of the system from four or three degrees of freedom to two degrees of 
freedom. In the present paper we treat for the first time an example, which only can 
be reduced to three degrees of freedom. We show to which extent our methods still 
work as before and how they reveal the dynamics of the deuterium in CDBrClF jl3j . 

For the already known wave function ( they are needed to construct the Hamilto- 
nian ) transformed from the number representation so as to be able to be plotted in the 
toroidal configuration space of action/angle variables jl4j we show that we can visually 
sort most of the states into ladders of states with similar topology. Each ladder has a 
relatively simple spectrum. Complexity arises from their interleaving. We also can rec- 
ognize underlying classical lower dimensional organizing structures by the fact that for 
these inherently complex functions the density ( magnitude squared ) is concentrated 
around them and the phase has simple behavior near them. The position of these 
structures in angle space reveals the nature of the resonance interaction causing the 
topology. It also allows the reconstruction, often by analytic use of the transformation 
back to generalized coordinates, which are near displacement coordinates, of the motion 
of the atoms which underlie this particular quantum state. By counting nodes in plots 
of the wave function density and phase advances in corresponding plots of the phase 
quantum excitation numbers can be obtained. In many cases, along with the polyad 
quantum number itself, we thereby obtain a complete set of quantum numbers for a 
state even though the corresponding classical motion is nonintegrable. These classifica- 
tion numbers can be interpreted as quasi conserved quantities for this particular state 
or for the ladder of states based on the same dynamic organizational element. Note 
there are usually several different types of organizing structures and hence ladders in 
the same energy region. 

Why are the wave functions that are so visually complex as to be unsortable in 
displacement coordinate space simple enough to be sorted in our semiclassical repre- 
sentation ? The answer is fivefold. First, when a spectroscopic Hamiltonian is used 
and as energy changes frequencies can move into resonance, that is rational ratio, a 
constant of the motion called the polyad number comes into existence ( see section 2 
). It can be used to reduce the dimension of the configuration space when the action 
angle variables, which come out in taking the semiclassical limit of the spectroscopic 
Hamiltonian, are used. Dimension reduction is always a great simplification. Second, 
the resonant form of the fitted spectroscopic Hamiltonian can be justified formally by 
using Van Vleck or canonical perturbation theory starting with a given potential sur- 
face |131ll5l[TB] in a way that takes in only the dominant dynamics. The terms that are 
left out for the most part are multiple small perturbations that in turn cause multiple 
small distortions to the wave function. This makes them harder to recognize but not 
fundamentally different. An excellent example of this is seen in P|, section 3 and quite 
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pictorially in |17j where figures 4 and 9 show the wave functions for the same states but 
computed from the full and the spectroscopic Hamiltonian respectively. The simplicity 
of the latter is enough to illustrate our point. The third reason for the relative wave 
function simplicity is that as the bound wave functions lie on a toroidal configuration 
space their parts must interfere constructively as each angle is changed by 27r. This 
coherence results in simpler but not always simple wave functions. The fourth reason 
is in the freedom to choose a particular type of reduced dimension canonical variables. 
These variables are designed with guidance from the polyad and the angle dependence 
of the resonance to have zero velocity when the phase locking caused by the resonance 
sets in. The wave function in reduced configuration space will then in turn further 
collapse onto the environs of ( in the sense of small oscillations about ) even lower 
dimensional surfaces or lines. We defer explaining and exhibiting this point till section 
3 which introduces needed notation. The fifth reason for simplicity is that, as alluded 
to above, the wave function is complex and can then be plotted separately in density 
and phase. Each part helps in visually sorting of wave function topologies and therefore 
dynamics and each yields the values of different quantum numbers. 

As organizational elements we do not necessarily use actual classical periodic orbits. 
In many cases it is simpler to study the general ( average ) classical flow and to relate 
to it some idealized guiding center and to use this idealized motion as the skeletal 
element on which the quantum density is concentrated. We can imagine that this 
idealized structure contains the essence of a large number of periodic orbits performing 
on average a similar motion. The quantum dynamics is not able to resolve all the 
fine details of the classical motion, nor the infinity of long periodic orbits, nor their 
enormous number of bifurcations. Its limited resolution only recognizes the average of 
the classical dynamics over cells of the size of its resolution. Therefore we interpret the 
quantum state as based on the motion of the atoms which corresponds to the idealized 
structural element. In many cases one state of the group of similarly organized states 
can be interpreted as the quantum ground state on this element and the other ones 
are various longitudinal and transverse excitations of the basic structure, hence the 
ladders. In some cases only sections of a ladder are found. 

Further in this paper, we see how the higher number and smaller spacing between 
observed states of a system with more degrees of freedom makes the whole investigation 
more complicated and in particular how the higher density causes more mixing of 
structures. This increase of structural mixing might be interpreted as one way in which 
effects of classical chaos enter the quantum behavior and results in a group of states 
being unassignable even in the semiclassical angle representation that brought 
simplicity in appearance to the density and phases of the assignable states. 

In section 2 we discuss the spectroscopic Hamiltonian, section 3 presents the semi- 
classical representation of the wave functions on the toroidal configuration space and 
section 4 contains the dimensional reduction. Section 5 explains the basic ideas how to 
extract the assignment of a wave function from its plot on the reduced configuration 
torus. Section 6 explains how from information about the dynamics in the reduced 
system we go back to reconstruct information about the original system. Section 7 
contains the classification of states and section 8 gives final remarks. 
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2 The model 



For the description of the chromophore dynamics of the molecule CDBrClF we use the 
model set up in ^HJ. It is based on 4 degrees of freedom as required by the nature of 
the observed overtone spectroscopy. They describe the three degrees of freedom of the 
D atom ( index s for the stretch and indices a and b for two bends, where index a is for 
the bend in the C-D-F plane and index b for the bend out of this plane ) and the C-F 
stretch ( index f ). The algebraic Hamiltonian fitted to experiment and to calculations 
on a fitted potential surface has a natural decomposition 

Hg = Ho + W (1) 

into a diagonal part 

-f^o = ^^jo]aj + ^ Xj^rna^ajolnam (2) 

j j<m 

(here the indices j, m run over the 4 degrees of freedom a, b, f, s) and an interaction 
part 

W = ksj^rniasOjal^ + alajajn]/(,2V2) + (3) 

'y[aaaaalal + alalabab]/2+ (4) 

SlaattaO^fa^f + alalafaf]/2 + (5) 

elababo'jra^j- + alalafaf]/2 (6) 

where aj and are the usual harmonic destruction and creation operators of degree 
of freedom j. In Eq.3 the indices j and m run over the 3 degrees of freedom a, b, f. 
For the coefficients we use the ones given in column 4 in table VIII in where we 
see that roughly ~ 2L0a ~ 2iOb ~ 2a;/ giving rise to the Fermi resonances anticipated 
by the terms of Eq.3 and the Darling-Dennison resonances anticipated by the terms of 
Eqs.4-6. Eq.2 has first an "oscillator" term and a second "anharmonic" Dunham term. 

The implied unperturbed basis of the space upon which Hg operates is the number 
basis 

|n >= \ns,nf,na,nb >= jn^ > |n/ > jn^ > \nb > (7) 

Tij are the quantum numbers or degree of excitation of the harmonic mode j and 
are eigenvalues of the number operator hj = ffljOj. \nj > is its pure oscillator mode 
eigenf unction. Hq is diagonal on this basis. The unperturbed ( first term of Hq ) 
energy of a basis function |n > is = J2j ''^j^j- When frequencies are near resonant 
as LOg ~ 2(jja ~ 2uib ^ 2ujf a polyad naturally arises with many consequences. First we 
can write = 2Pujb with the polyad P defined as 

P = ns + {nf + na + nb)l2 (8) 

with the corresponding operator given as P by replacing rij — > fij in Eq.8. Clearly all 
basis functions with rij giving the same value of P are degenerate. If basis functions 
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are sorted by polyad, the Hamiltonian matrix which can be evaluated algebraically, is 
expected by perturbation theory ( degenerate basis functions interact most strongly 
) to be nearly block diagonal suggesting that P is a near constant of motion. In 
fact, it is totally block diagonal as P is here an exact constant of the motion since P 
commutes with Hq by construction albeit not with the true Hamiltonian to which 
Hq is, as said, to a high order of perturbation an approximation. When each block is 
diagonalized, each eigenfunction 



in block P has a good conserved quantum number P which roughly measures via the 
equation for £"2 the degree of excitation of the state in units of the lowest frequency. 
One can now study each polyad separately ( second consequence of P ). Here we will 
concentrate on the one 161 state polyad P = 5 which is the most complex one in the 
range of measured energy. 

3 Semiclassical representation 

Having the {c^} now allows us to immediately obtain, without any extra work, semi- 
classical eigenfunctions in action/angle variables ( good to the lowest two orders in h 
). We do this because as we will see such eigenfunctions can be viewed due to the exis- 
tence of the polyad as effectively having one less dimension. Hence the problem must 
now be transformed from the number representation to a semiclassical action/angle 
representation. The drop in dimension is characteristic of classical mechanics where a 
constant of the motion, here the classical analogue of P, can by canonical transforma- 
tion allow one to reduce the dimension of the phase space. Here the reduction is from a 
four degree of freedom system with action/angle coordinates /, to a reduced system 
with specific value of P and three degrees of freedom with action/angle coordinates J, 
ijj. This of course will accomplish the first promised simplification, that is dimension 
reduction. 

To make a transition to the corresponding classical Hamiltonian we bring the oper- 
ators in Hq first into symmetrical order, and then apply the semiclassical substitution 
rules 



where Ij and (pj are the classical action and angle variables. The resulting classical 
Hamiltonian is using Eqs.2-6 




(9) 



nGP 




(10) 



^0 = - V2)c^j + E - V2)(/m - 1/2) 




j j<m 



(12) 




(13) 
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6IaIfCOs{2^a-2(l)f) + 

eIbIfCOs{24)b - 2(pf) 

It has the conserved quantity 



(14) 
(15) 



K = Is + [Ia + h + If]/2 (16) 

Because of the symmetric ordering required for the transition to classical dynamics the 
zero point between the eigenvalues of P and the corresponding value of K is shifted by 
5/4. 

The reader, as an aside, might wish to contemplate how easily this transition was 
taken. For the more familiar H = T+V Hamiltonian form one must solve the Hamilton- 
Jacobi equation to get such variables; indeed a difficult process even if the system is 
completely integrable. The simplicity is due to the fact that, as said, this should be 
rigorously done starting with a potential surface and using perturbation theory. A price 
is paid for circumventing perturbation theory and using a fitted spectroscopic Hamil- 
tonian. This price is that the relation of the and aj to the generalized momentum 
and position variables pj and qj is rigorously unknown. We will return to this issue 
later when going back from reduced space to displacement coordinates; our so called 
lift process. 

If one now quantizes any single harmonic oscillator term in Eqs.11-15 using Schroedinger 
correspondence Ij — —id/d(f)j, (f)j — (f)j it is clear that, since the angles define a now 
toroidal configuration space, 

\nj) — exp{inj(f)j), |n) — exp(m • (f)) 

with Tij corresponding to Ij — 1/2. From Eq.9 we get 

*f ^ *f (^) = E 4 exp(m • 
HeP 

That is, the semiclassical wave functions are already known and given by Eq.l8. 

4 Dimension reduction 

Dimension reduction can now be carried out by noting that here, using the polyad of 
Eq.8 to eliminate the fast degree of freedom s, that the j-th eigenstate can be written 

^f{$) = eMiP<t^s) E 4',n.,n„n, H exp(mfe(<^fe - <^,/2)) (19) 

na,nb,nf k=a,b,f 

= exp(iP(/),) J2 '^P,na,n,,nf li exp(mfcVfc) (20) 

na,nb,nf k=a,b,f 
= exp{iP(f)s)Xj {ipa, tpb, tpf) (21) 



(17) 



(18) 
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This shows that the eigenfunctions in a given polyad P have a common phase factor 
dependence on 0^ and really only depend on which is a function of the three angles 

il^j = - ct>s/2 (22) 

for j = a, 6, /. These three angle make up a three dimensional toroidal configuration 
space upon which Xk is situated. In fact Xk J^^^ Fourier decomposition of the 
reduced dimension wave function on the configuration torus. More rigorously the wave 
function xj" defined by Eqs.l9- 21 actually is the semiclassical wave function after the 
classical constant of motion K is used to reduce dimension by a canonical transforma- 
tion. Using the angle transformation of Eq.22 suggested above by the existence of P 
plus the trivial transformation 

e = ^s (23) 
the canonical transformation permits that the new actions are K, given by Eq.l6 and 

J, = (24) 

again for j = a,b,f. 9 is the cyclic angle. Staring with equations 11 to 15 and 
using equations 16, 11-24 a new Hamiltonian in the six canonical variables J, ip albeit 
parametric in P or can now be written so directly that to save space we omit it. 
The torus is the configuration space. Using the Schrocdingcr quantization for the 
new variables is only semiclassical, i.e. is correct in the two lowest orders in h. Its use 
gives the semiclassical expression 

Xj {'4^aAbAf)= X] ^P,nf,na,nb exp(i(naVa + WfeVb + "/V'/)) (25) 

for the eigenfunctions on the reduced configuration space. 

In the new coordinates the basis functions with polyad number P and three indices 
nf,na,nb corresponds to the product of three plane waves or rotors as 

Yl exp(infeV'fe) (26) 

k=a,b,f 

Each single plane wave factor loops the in the positive direction. The rotor is that 
into which the single mode harmonic oscillator has been transformed. Since Eq. 25 is 
a manifestly complex function it can be written as 

= \^^^f{^|;)\e'^f(^^ (27) 

where the phase function is given as 

$f = arctan(/m^'f /Ee*f ) (28) 

For the basis function of Eq.26 it is well to note for future reference that the density is 
constant and the phase is 

= rifi/jf + ria'^a + ribtpb (29) 
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Symmetry properties of the system in the toroidal configuration space will be used in 
the following. The Hamiltonian is invariant under a simultaneous shift of all angles ipj 
by TT. Therefore the reduced configuration torus covers the original configuration 
space twice and all structures show up in double, even though they really exist only 
once. In addition the system is invariant under a shift of any angle ipj by 27r and it is 
invariant under a reflection where all angles go over into their negatives. 

5 Methods to evaluate wave function plots 

It is seen that the simple basis functions have constant density and any eigenstate dom- 
inated by a single basis function has a density without sharp localization. Resonances 
as for example 2L0a ~ ujs will be seen to cause localization about a line ipa = constant. 
This follows as ipa = 4'a ~ = constant when differentiated with respect to time 

gives the frequency relation. Hence it is seen that resonances are associated with local- 
ization and the fact that dtpj/dt = 0. It can be claimed that by using angle coordinates 
that slow to zero velocity at resonance we here assure that the wave function will "col- 
lapse" onto and about a lower dimensional subspace called the organization center. 
This in reverse gives a way to recognize the influence of resonances, namely localization 
of the wave function on the configuration torus. Each linearly independent locking of 
angles therefore reduces the dimension of the subset of configuration space in which 
the wave function is concentrated by 1 dimension. 

The fact that the wave function is on T^, a three dimensional torus, presents graph- 
ical and visual challenges when attempting to sort wave functions by inspection. To 
meet these challenges we employ the fact that is 1:1 with a cube with sides of angle 
range 2tt and identified opposite boundary points. That is a point on one surface of 
the cube is the same as a point on the direct opposite side of the cube. The density 
and phase are plotted in such cubes. The density plots now over the interior of the 
cube are organized over the whole cube or about interior planes or lines or points if 
zero, one, two or three independent resonant couplings respectively are active. Over 
(or "along" in the case of a line) the organizing structure a smooth density with no 
sharp localization should appear. 

Nodal structures will be visible and countable in directions perpendicular to the 
organization element and clearly will be associated with a localized direction. The 
count of such transversal nodes supplies for each direction of localization a transverse 
quantum number t, that replaces an original mode quantum number n which has been 
destroyed by the resonant interaction. The wave function in all the localized directions 
can be considered as qualitatively similar ( a continuous deformation ) of an oscillator 
state of the corresponding dimension. By the superposition of the running wave basis 
functions with appropriate weights and signs locally around the organization center the 
nodal pattern of an oscillator is reproduced. The transverse quantum numbers t^ are 
the corresponding oscillator excitation numbers. 

The phase functions <I> according to Eq.28 are smooth and close to a plane wave 
in the neighborhood of the organization center. They tend to have jumps by tt and 



8 



singularities away from the center. Therefore phase advances along fundamental cycles 
( loops on the torus or lines that connect opposite boundary points on the cube ) 
of the organization center are well defined and are necessarily integer multiples of 2tt. 
They provide longitudinal quantum numbers 1^ which also replace the interacting mode 
quantum numbers. 

Some states present what at first look appears to be a familiar but unexpected phase 
pattern in that the phase picture is too simple. Cases will appear when the density is 
localized but the phase picture has the globally striped flag appearance of the nearly 
noninteracting mode case. The strangeness disappears if ideas from nonlinear dynamics 
are employed. Nonlinear dynamics tells us that at low P most of the motion takes place 
on what is called the primary zone of phase space. And the now interacting states 
are quantization of tori in this region. This zone has no central organizing structure 
other than the whole configuration space itself. It does contain quasiperiodic orbits 
corresponding to the uncoupled modes. 

As P changes coupling of a subset or even all modes can occur by at least two 
mechanisms. First is one where the existing mode periodic orbits bifurcate and create 
new secondary zones based with newly created organizing structures ( periodic orbits, 2 
tori, 3 tori e.t.c. ). Also new secondary zones may be created seemingly spontaneously 
due to what are called saddle-center bifurcations. In these secondary cases the motions 
will have changed and states quantized based on tori or near tori ( think about the 
resolution of ^ ) will have density and phase diagrams as described above with (P, tk,lk) 
quantum numbers. 

The second mechanism of coupling which will manifest itself in density localization 

is by continuous deformation; a process visualized by mentally drawing the primary 
configuration on rubber sheets and stretching and pulling the sheets in various direc- 
tions till the localized form is achieved. Under such continuous deformation the phase 
function $ of the continuous, single- valued complex wave function *(V') always must 
change by an integer multiple of 27r as its argument goes around any closed loop 7 in 
configuration space. In our case prior to deformation ^ = n ■ ^ and that integer is Uj 
if the loop is in ipj ( this loop will be called 7^ ) i.e. as ipj — >■ ipj + 27r, + njlir. 

More generally for any $ this can be stated as the loop integral relation 



where rrij is an integer and A is the deformation parameter. Now under a continuous 
deformation it is clear that rrij or in our case rij must remain constant as a continuous 
integer valued function. In our case this rij can be read off from the phase diagram as 
equal to the number of revolutions of the phase along the loop 7^- achieved as ipj 
tpj + 27r. For the noninteracting modes the afore mentioned quantum number I is just 
the n. 

Prom the point of view of basis state mixing even though the phase is behaving as 
if a single basis function is totally dominant, significant state mixing can occur. This 
can and will occur in proportion to the amount of localization of the density. Weak 




(30) 
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localization tends to occur when one basis function is dominant, the one with the same 
n. 

The fi assignment holds till a small change in A causes an abrupt change in 
that is when singularities occur. At this point continuous deformation will stop. Now 
since the wave function must be continuous, singularities in phase can only occur where 
the density and hence the wave function is zero. For us this occurs at the nodes and 
perhaps the outer edges of the localized oscillator. In short when nodes occur one can 
no longer use this simple full set of n quantum mimbcrs. One can still count phase 
revolutions along loops inside an organized region ( parallel to the organizing center 
) to get the I quantum numbers. Great care should be exercised in deciding if for a 
localized density the phase is a continuation of a plane wave or not. To do this ^{tp) 
must be a smooth function in the neighborhood of the loop along which one is counting 
the phase revolutions. 

6 Ladders and motions in displacement coordinates 

Once the wave functions in their semiclassical reduced dimension representation are 
diagrammatically presented and quantum numbers assigned as described in the previous 
section, as will be seen in the next section, most often states can be visually sorted 
and ordered into ladders of states with a given topology. Each ladder tends to lie in 
the Husimi sense in one of the phase space resonance zones discussed previously. The 
simplicity of sorting is due to the fact that as said an idealized classical organizing 
element underlies each ladder. This structure can easily be written as a configuration 
space relation between the tp variables. The structure will allow us to immediately 
determine what type of resonance interaction is creating this structure. An approximate 
procedure that transforms the variables ^ back to displacement space reveals the atomic 
motion ( modes of motion ) whose quantization yields the ladders of states. This process 
we call "lifting" and refer to it as the "lift" . We now describe the lift in a way that starts 
with the least approximate and most difficult to implement algorithm and proceeds to 
simpler to implement ones that still retain the qualitative spirit of the motions. 

To determine the key interactions we only illustrate the procedure with a simple 
example. Assume that the organization structure is a line ( or fiber ) along ipa = 0. 
Since ipa = fpa — fps/^ then ipa — 4>s/^ ~ "^a ~ i^s/^ ~ and the 2cUa = phase locking 
relation speaks to us of a 2:1 Fermi resonance between the bending mode a and the 
stretch mode s. Other cases are similarly treated. 

Now for interpretation of the motion a trajectory in the (J, V") phase space need be 
found which when projected into configuration space should represent the organizing 
element. In reverse the knowledge of the organizing element is helpful in finding an 
appropriate trajectory. To start this search estimates of the actions will be useful. For 
uncoupled modes the assigned nj quantum numbers using Ij = rij + 1/2 suffices. To 
obtain estimates of the initial actions the quantum mechanical average of the Jj can 
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be used and trivially computed from the known wave functions as 



< Jj >=< *l - iwrl^ >= E 




(31) 



The ladder state that is most localized should be used for this and to check an 

estimate the variance can be used. With this wc now have initial conditions (< Jj > 
and a value of ^lJj on the organizing structure) with which to uses Hamilton's equations 
to get tpj{t), Jj{t). Then first we reconstruct the cyclic angle 9 by integrating 



Next we undo the canonical transformation to get the corresponding (f>j{t) and Ij{t). 
Then wc assume some idealized harmonic model for the elementary degrees of freedom 
and form the displacement qj{t) and its conjugate momentum (i) as 



The functions qj{t) will be interpreted as representatives of the atomic motion. 

Diagrams in which the qj are taken as coordinates with t as parameter give the 
more conventional picture of the atomic motions. In the simplified methods the av- 
erage values of Ij which giving the scale of the qj motion will replace I{t) in Eq.33. 
These methods arc simple because again Eq.33 is used but the time dependence and 
therefore the computation is suppressed and the relations gotten between the qj from 
the organizing structure equation that relates cpj along with the use of Eqs.23, 24. 

In preparation for section 7 let us discuss a few very simple cases. First there is the 
possibility that the density is distributed rather homogeneously over the whole and 
the phase function comes close to a global plane wave with 3 component wave vector 
{na,rH),nf). This means no locking at all, all modes run freely with their individual 
frequency and the excitation numbers of the various modes are rij read off from the 
wave vector for j = a,b, f (i.e. count the revolutions as one traverses the phase diagram 
in the j'th direction) and Ug = P — [na + 725 + nj]/2. If the density is concentrated 
in a plane, e.g. the plane V'o = constant then the angle Va does not move in the 
long run. It only may oscillate around the value of the constant. ^Prom the canonical 
transformation of Eq.22 we see that this means 0a — 0s /2 = constant for the angles 
of the original degrees of freedom. Therefore we have frequency and phase locking 
between the modes a and s. The other modes b and f run freely and their excitation 
is according to the wave vector of the phase function in the plane of high density. 
Analogous considerations hold for other planes of concentration of the density. If the 
density is concentrated along a line, e.g. the line ipa = constant and = constant 
( this is a fiber running in f direction ) then we see from Eq.22 that 0a — 0s/2 and 
06 — 0s /2 are kept close to constant values and this means locking between a and s 
and simultaneously locking between b and s. Of course, then also a and b arc locked 
relatively to each other. The mode f runs freely in this case and its excitation number 
follows from the wave number of the phase function along the fiber of high density. 
Analogous considerations will be applied for a concentration of the density along other 
1 dimensional lines in reduced configuration space. 




(32) 



Q^i) = \/2^cos(0,(t)),p,(t) = y2i^(i)sin(0,(t)) 



(33) 
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7 The motion behind individual quantum states 



For polyad 5 with 161 states, the given transformations and equations lead to a three 
dimensional toroidal configuration space which here is a cube with coordinates tpa, tpb 
and il^f. Three dimensional plots of the semiclassical state functions look too often like 
large globes of ink and were not informative. Here we will resort to cuts of the cube ( 
e.g. cut tpa = V'b) plot density and phase in such cuts. 

The question is now where to cut. Clearly a knowledge of the classical flow would 
be most helpful here. In it's absence analysis as follows is also useful. The cuts should 
contain organizing elements about which the density is localized. The phase should be 
simple over the regions of the plane having high density. Note the idea of an organizing 
structure is that one has found the confi,guration part of the quasiregular resonant zone 
of the coupled modes. The phase simplicity can be used as a tool for finding organizing 
structures if they exist. Here we have used both trial and error as well as dynamics 
guidance to find the cuts. 

Dynamics guidance means for example that in CDBrClF the resonance with closest 
frequency is LOa ~ and its associated ka^b in the Hamiltonian is largest. This sug- 
gests that modes a and b should be coupled in the plane V'a = V'6 + constant which 
assures uj^-^^ = '^^^ ■ Note the constant will vary as explained below in a dynamically 
meaningful way, i.e. it will be used for state classification on a ladder of states. 

The Fermi choice a;^ = 2iOa along with oog = 2a;(, would then give by similar reasoning 
the organizing line ipa = ci, Vfe = C2. For some states several organizing structures 
could be used and the same dynamics revealed. Our experience is to choose those 
corresponding resonances which seem more important in H and which give longer 
ladders of states. 

7.1 Scheme for a leirge number of states throughout the polyad 

Starting at the bottom of the polyad we find about 64 states that definitely lie in the 
Darling Dennison coa = U)b resonance class in that they have densities localized about a 
plane ipa = ''Pb + constant albeit often with different constants. They also have simple 
phase plots in these planes. At the lower end a few of them could also be organized 
by the a;^ = 2a;a, Ug = 2oJi, Fermi resonances , i.e. about a line in f direction; but they 
occur when stretch excitation is zero and this classification is not very physical. 

Some of the lowest states could also be assigned for all modes by quantum mimbers 
nj and correspond to continuous deformations of basis functions. Although these nj are 
useful in the lift for obtaining actions these assignment misses the phase and frequency 
locking and the localization of the density. In addition we looked very carefully for 
primary tori in the classical phase space at the corresponding energy and did not find 
any. Therefore EBK quantization cannot be applied. And classical dynamics shows 
that the assignment by nj for j = a,b, f is unphysical. 

To illustrate our ideas let us consider two states, state s = 4 and s = 46 at the 
bottom and the middle of the ladder respectively. The former state is localized about 
the plane ipa = V'ft the latter about V'a = V'6 + ""/^ as seen in figures la and Ic for 
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the former and figures 2a and 2c for the latter. Both locahzations give the DD coa = tOb 
expectation. Note how only in the former case also the line V'a = 7i"/2, V'6 = 7r/2 
could have been used as organizing center. In the diagrams all features come in double 
because of the symmetry mentioned at the end of section 4. Figures 1 and 2 show that 
if one stays on the organizing plane and heads say in the f direction the torus is looped 
during a 27r change of ipf- This says that ipf motion is free and that on the phase 
plot following any path that along which -0/ — > ^/ + 27r the free exp(m/'0/) factor in 
the wave function should go to exp(in/(V'/ + 27r)). In Fig. lb or Fig. 2b respectively 
asV/~^V'/ + 2vr and if we count the revolutions along a line in f direction one sees 
a phase advance of 2-k and IOtt respectively. We conclude that n/ = 1 for state 4 ( 
Fig. lb) and n/ = 5 for state 46 ( Fig.2b). Sometimes here n/ will be usefully called If 
to emphasize that it is a longitudinal quantum number obtained by following the path 
of highest density around the f loop. 

The phase change accompanying a loop of the a direction in figures lb and 2b by 
the nature of the cut causes a looping on the V'o = "06 + constant plane in the complete 
configuration space. This looping advances both and iph by 27r and on this plane 
the wave function must cliangc as ccxp(m • tp) — > X] cexp(i(n'i/; + 27r(na + nj,))- The 
choice of the ipa = "^Pb^ constant plane required that we see a simple plane wave factor 
exp{i{natpa + ribipb)) which means that all the important basis states with large mixing 
coefficients must have this type factor. Hence the wave function phase must advance 
as we sweep ipa in this organizing plane by 27r(na + Ub). We denote this by la+b, I for 
longitudinal, i.e. motion along the organizing center as opposed to transverse to it. 
This longitudinal quantum number clearly tells us the total number of quanta in the 
Darling Dennison lock. 

For state 4 in figure lb la+b = 9 and for state 46 in figure 2b la+b = 5. Using 
the polyad expression for the stretch mode we get ng = P — {nj — la+b)/'^ = 0- Now 
a partial ( three out of four needed quantum numbers ) assignment can be made as 
{P-,lf = nf,ns) or equivalently {la+b, If = nf^ng). State 4 is now (5,1,0) or (9,1,0) 
respectively and state 46 is (5, 5, 0) in both schemes as recorded in table I. In principle 
na and no longer exist as good quantum numbers except when as explained these 
numbers represent conserved integer loop integrals as in the continuous deformation 
case for state 4 where both phase diagrams ( Figs lb and Id ) look as they represent 
2D plane waves. Figure Id then tells us that = and 725 = 9. This is useful only for 
obtaining /a = na + 1/2 = 1/2 and lb = nb + 1/2 = 9/2 which will in turn be useful in 
the lift process. 

For state 46 phase simplicity only happens in the organizing plane. Fig. 2d shows 
no useful phase information. The plot shows jumps along lines and ramp singularities 
where the phase value depends on the direction of approach. 

At this point and in a similar way a full 64 states for P = 5 can be assigned by ng and 
nf. Actually another 32 seem to resemble this picture. The resemblance is only seen 
as trends from the systematic following of the 64 other ones. For these increasing but 
still smaller effects of many other resonances are taking their toll and they significantly 
degrade the assignable states into unassignable more heavily mixed quantum analogues 
of classical chaos. The table I indicates which states have less certain assignments. 
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Certain of the states will be marked in the table as demixed states and given in 
s/s' notation, i.e. 50/51. For such pairs it will be seen that the eigenstates involved 
are nearly degenerate and usually do not have the features of any of the states on any 
of the ladders we identify. The near degeneracy is an indication on that possibly two 
( noneigen ) states exist with even more degenerate energies that do fit on ladders 
but fail to appear due to mixing by accidental degeneracy. By assuming for such 
eigenstates various weighted mixings the demixed dynamically characterizable states 
can be recognized and included in the ladders. 

For further simplifying and dynamical classification the states of the DD ladder can 
now be sorted by common values of rij and Ug ( or equivalently la+b and // ), see table I. 
Clearly one quantum number is missing. In a lower dimensional problem ( e.g. ^) in a 
similar case of DD lock ( there between local stretch modes n and m ) the longitudinal 
In+m which was determined by counting phase advances along an organizing line in the 
V'm = V'n direction, was complemented by counting nodes in the direction transverse 
to the organizing line by using this number for a transverse quantum number t. Here 
only remnants if any of these nodes are seen. The greater mixing in higher dimension 
wipes out this clarifying element. What does remain common is that states with the 
same other quantum numbers as above can be ordered by how far significant density 
extends into the direction transverse to the organizing plane. For the DCO this gives 
the same t ordering. 

Here we adopt this criterion to organize the states in the //, la+b or nj Us sub- 
ladders and also observe in Table I that this ordering simultaneously gives an energy 
ordering. This t ordering number is our final quantum number. Additional support for 
our classification comes from the classical mechanics where trajectories at eigenstate 
energies flowing along the organizing plane fluctuate further in the transverse direction 
the greater is t. Most gratifyingly each subladder has an energy spectrum that is quite 
reminiscent of an anharmonic oscillator. As such locally and physically for each sub- 
ladder of fixed // and la+b but varying t we envision a dynamic potential causing the 
system to remain near the organizing structure. The levels of this potential are the 
states of the subladders. 

Even at this stage where other ladders discussed below are omitted spectral com- 
plexity arises from the interlacing of energy levels of different subladders. The roughly 
96 states on the ladder which only demands the DD-a/b lock for its construction gives 
a complex spectrum and nontransparent energy spacings. 

The reader will note that some ladder states in table I are missing. These states 
could not be found but their disappearance is confirmed using the anharmonic models 
expected energy spacing. Formally these states lie among the highly mixed ( chaotic ? 
) states we discuss below. They are mixings of more than two dynamically identifiable 
states and do not show the near degeneracy that our demixing process identifies. 

The actions Ij can be gotten for the modes s and f using Ij = nj + 1/2, in addition 
we have la + h = la+b + 1- Then = J a and = Jb is gotten from the state functions 
using Eq.31. To observe approximate averaged motion of the Deuterium atom the Ij 
can be used in Eq.32 for 9 to enable us to get I{t) and (j){t) and then using Eq. 33 to 
draw qualitatively pictures for the various pairs {qj,Pj)- Actually a great deal can be 



14 



said without the integration for 9. In (g^, qa/h/f) planes a roughly rectangular region is 
seen whose sides are proportional to \/2l7 and The deuterium trajectories 

motion interior to the rectangle will be quasiperiodic motion. How the exact trajectory 
depends on the initial choices of 4>s and '^a/b/f 1^ really not important. In the 
{qa,Qb) plane organization structure leads to ellipses whose relative ranges ( excursion 
from zero displacement ) in the {qa, qb) variables is \/Ia/Ib- The ellipses eccentricity 
depends on the organizing structures phase shift, viz. = fpb + ^ ^ <t>a = <t>b + ^■ 
(5 = and tt give zero angular momentum and something approaching a straight line 
motion while 5 = ib7r/2 gives maximal angular momentum and elliptical motion which 
becomes circular if /„ = I^. The t value is the out of phase motion and causes the 
trajectory to change slowly its eccentricity. 

7.2 Scheme for the upper end of the polyad 

At the upper end of the polyad we take to start states 160 and 161. A cut in = ( 
or any other value of ■06 ) reveals as we see in Fig. 3c for state 160 what looks like two 
columns of density localized around = vr and -i/'/ = tt. In comparison state 161 ( not 
shown ) only has one column. A cut in the plane ^6 = constant is required to see that 
the columns really exist. Figure 3c shows that they do for state 160 and it also holds 
for state 161. Clearly for both these states the line = tt, ■0/ = tt is the organizing 
structure. Since it loops in -06 mode b, is decoupled and the phase picture for both 
states shows ( figure 3b for state 160, state 161 not shown ) that as -06 — > V'fc + ^vr along 
the organizing line no wave fronts are crossed and hence there is no phase advance. As 
such lb = nb = 0. For the upper end of the polyad also classical dynamics shows a flow 
in b direction. 

Now the ipa = "ipf = organizing structure implies that the phase locks are 0a = 
0s/2 + TT and 0/ = 0s/2 + 7r. Time differentiation of these relations show that the 
frequency locks are of 2:1 Fermi type with the stretch mode locking with the bend 
mode a and the mode f. Clearly for the organizing structure ^pa = ^/ or equivalently 
then 0a = 0/ is a valid phase lock and Wa = is a DD 1:1 frequency lock. Since 
all three types of terms appear in the Hamiltonian we can be assured that all three 
frequencies are in lock and that the quantum numbers na, nf and no longer exist. 
They must be replaced by three quantum numbers of the lock one of which can be P 
and the other two can be taken as the number of nodes seen in the V'b = cut along 
the antidiagonal as = 1 and along the diagonal as t2 = 0. For state 161 since one 
column exists ti = ^2 = 0. 

Table II shows the assignment for the states whose analysis is based on the same 
organizing structure. Figure 3a shows a vertical cut into the cube along the antidiagonal 
of figure 3c. It reveals at -00 = tt a vertical nodal line corresponding to the node on 
the antidiagonal line -0/ = —ipa in figure 3c. The assignment of states 161 and 160 in 
table II is complete. P — nb is the number of excitations in the three way lock and 
the quantum numbers ti and t2 tell us how much of this excitation is transverse to the 
columns and therefore is a measure of the extension of the density outward from the 
organizing line created by any two of the three resonances mentioned above, only two 
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of them are independent. 

All entries in table II use the same organizing structure as these states. Two of them 
needed to be demixed before being assigned. They are not continuous deformations 
of basis states ( the phase functions never look like plane waves ), only the action 
can be obtained as /(, = n(, + 1/2. Estimates of la, If and Is must be gotten from 
computing the quantum state average of Ja, Jf and by the canonical transformations. 

An idea of the internal motion of the deuterium and the CF stretch can be obtained 
from the actions, phase relations and the ti and t2 values. In all {qi,qj) planes the 
actions again define the maximal displacement ^/2Tj of each qj from its equilibrium. 
Clearly this range is meant in the sense that an average is made over the values of the 
two variables not represented in the plane. The planes {qb,qj), j = o-jfjS will show 
rectangularly bound quasiperiodic motion. In the {qs,qa/f) plane a U shaped region is 
swept out. Here q^/f reaches its extreme values as qs reaches its maximum displacement 
and qs sweeps its range twice for each sweep of qa/f- Increasing ti and/or t2 causes the 
lifted trajectory to oscillate in a tube about this basic U lift. Of course qf and qa reach 
their maxima ( minima) in phase as the three modes s, a and f are phase locked. 

7.3 Patterns in the dense region of the polyad 

In the middle of the polyad there are approximately 50 states with very complicated 
wave functions which we could not classify. It is possible that further and closer analysis 
could reveal some systematics. These states are dispersed among the states mentioned 
in subsection 7.1. But interestingly a few simple states which do not fall into any 
scheme recognized are dispersed in this region. In this subsection we present a few 
such states. 

State 111 has its density concentrated in the plane ipf = ipf, in which plane the 
phase function also comes close to a plane wave with la = and = 10 ( see parts 
a and b of figure 4 ); this simple n,, = 0. In the density plot we see a strong tendency 
towards a-b coupling but it is not perfect. The parts c and d of the figure show density 
and phase in the transverse plane ipf = tt/2. We see a very good concentration of the 
density along the line ipb = ipf indicating t = with respect to this organizing plane. 
Therefore in this state all excitation is in the in phase combined DD-b/f ( bend of 
deuterium and stretch of fluorine ) motion. The degrees of freedom a and s are not 
excited and the assignment is complete. 

Numerical results for state 103 are shown in figure 5, part a and b show density and 
phase in the plane ipf = which is the plane of high density and thereby the organizing 
structure. Part c and d show density and phase in the transverse plane V'a = tt/S. We 
clearly see how the density is concentrated around the planes V"/ = or tt. In the plane 
of high density the phase function comes close to a plane wave with = 1 and = 4. 
A few other states ( e.g. 108 and 148 ) fall into the same scheme. For the lifted motion 
in the displacement space we can draw the following conclusions: a and b both run 
freely with an excitation corresponding to their loop numbers read off from the phase 
function in the plane ijjf = 0, i is completely coupled to s in a 2:1 Fermi resonance. 
P=5 and the t = from figure 5c complete the assignment. 
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State 140, Figure 6, is an example with a concentration of the density around a 
point, in this particular case the point ipa = — vr/2, ipb = vr/2, il^f = 7r/2. This means 
coupling of all degrees of freedom at a particular phase value. The phase function 
is rather simple globally and indicates n,j = 0, Ufy = 0, Uf = 8. Because of polyad 
conservation we then conclude that = 1. This means that all four actions are 
Ij = Hj + 1/2. Accordingly we have a coupled 2:1 Fermi-s/f motion whereas the formal 
coupling of a and b with s is rather irrelevant since their excitation numbers are 0. This 
interpretation is based on the phase function even though the density is concentrated 
around a point, where logically we should classify states by transverse excitations in 3 
directions ( are all zero in this particular state). 

8 Conclusions 

We have shown how our previously developed methods can also be applied to systems 
with 3 degrees of freedom after reduction. For seventy percent of the states the strat- 
egy works successfully. We investigate the wave functions constructed on a toroidal 
configuration space. We investigate on which subsets the density is concentrated and 
evaluate the phase function on this subset of high density. This allows us to obtain a 
set of excitation numbers for the states; in the best cases we obtain a complete set of 
excitation numbers which provide a complete assignment of the state. The excitation 
numbers can be interpreted as quasi conserved quantities for this particular state. 

The success of the method is based on 3 properties of our strategy which more tra- 
ditional approaches starting from potential surfaces do not have. First, we start from 
an algebraic Hamiltonian which naturally decomposes into an unperturbed Hq and a 
perturbation W. Perturbation theory assures us that W contains only the dominant 
interactions, leaving out those that do not affect the qualitative organization of the 
states. In addition having the ideal Hq results in the classical Hamiltonian being for- 
mulated from the start in action and angle variables belonging to Hq. For a potential 
surface it is very complicated to construct an appropriate Hq and the corresponding 
action/angle variables. Our knowledge of the correct Hq allows the continuous defor- 
mation connection of some cigcnstates of H with the corresponding eigenstates of Hq 
which was useful for many states. Second in the algebraic Hamiltonian the polyad ( 
which for the real molecule is only an approximate conserved quantity ) is an exactly 
conserved quantity and allows the reduction of the number of degrees of freedom. We 
always deal with one degree of freedom less than the number of degrees of freedom on 
the potential surface. Third, our strategy depends to a large extent on having a com- 
pact configuration space with nontrivial homology where the phase of the wave function 
is important and where the phase advances divided by 2tt along the fundamental cycles 
of the configuration space provide a large part of the excitation numbers. In part also 
feature 1 is responsible for the importance of the phases. This is an essential difference 
to working in usual position space. In total we have found a semiclassical representa- 
tion of the states which is actually much simpler to view than the usual position space 
one. 
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The classical dynamics of the system shows rather little regular motion and is 
strongly unstable. Then the question arises in which way the quantum states reflect 
anything from the classical chaos. We propose the following answer: In the middle of the 
polyad a large number of states do not show any simple organization pattern. For some 
of those states it might be that we simply have not found the right representation in 
which to see the organizational center. However we doubt that this is the case for many 
states. For many of the complicated states the wave function is really without a simple 
structure, it seems to be a complicated interference between several organizational 
structures. Then immediately comes the idea to mix these states with their neighboring 
states to extract simple structures in the way as we did it for states 50 and 51. For most 
of the complicated states this method does not lead to any success. This indicates that 
underlying simple patterns are distributed over many states so that demixing becomes 
rather hopeless and maybe also meaningless. In the spirit of reference |L9j this indicates 
that the density of states is already so high that many states fall into the energy interval 
over which states are considered almost degenerate and mix strongly. This is a strong 
difference to the systems we have investigated before which could be reduced to 2 
degrees of freedom. There the density of states was much lower such that mixing only 
occurred occasionally between 2 neighboring states. This stronger mixing gives some 
clue as to how the quantum wave functions become complicated for classically chaotic 
systems and for sufficiently high quantum excitation. 

Finally there may come up the question why we do not use some of the common 
statistical procedures ( like e.g. random matrix theory ) to treat our system when 
effects of classical chaos become evident in a quantum system. The answer is: We 
are not satisfied with only generic (universal) properties of the system or with the 
investigation to which extent such universal properties are realized by our system. We 
aim to investigate each individual quantum state one after the other and to work out its 
individual properties in order to arrive in the best case at a complete set of quantum 
numbers in addition at a detailed description of the atomic motion underlying the 
individual state. Previous examples have shown that our procedure provides a complete 
classification for simpler systems which can be reduced to 2 degrees of freedom. As the 
example presented in this paper shows, this program is doable with good success for a 
large number of states even in a system with originally 4 degrees of freedom, reducible 
only to 3 degrees of freedom and with a complicated classical dynamics. 
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Table I 



Table I: This table gives the classification and assignment of all states organized around 
planes i^b = i^a + constant. The first column gives the number of the state ordered 
according to energy. If the structure only becomes clear after demixing, then it contains 
in addition the number (s) of the state (s) with which it is demixed. The second column 
gives the energy in units of cm~^. The third and fourth columns give the two longi- 
tudinal quantum numbers If and la+b- The fifth column gives the transverse quantum 
number t. For such states where the individual quantum numbers na and ra^ exist, 
these are given in columns 6 and 7. If these phase functions show some first indications 
of developing singularities, these quantum numbers are put in brackets. States where 
the pattern is not so clear but where the general trends suggest the assignment have a 
* in the last column. The total assignment is (P, or n/, or n^, t). 



Number (s) 


Energy 




la+b/ris 


t 




nb 




1 


8908 





10/0 








10 




2 


8997 





10/0 


1 


1 


9 




3 


9080 





10/0 


2 


2 


8 




5 


9156 





10/0 


3 


3 


7 




7 


9226 





10/0 


4 


4 


6 




10 


9290 





10/0 


5 


(5) 


(5) 




13 


9348 





10/0 


6 


(6) 


(4) 




16 


9400 





10/0 


7 








20 


9446 





10/0 


8 








25/24 


9485 





10/0 


9 








29/28/30 


9526 





10/0 


10 






* 



to be contiunue on next page 
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Number(s) 


Energy 


If = Uf 




t 


na 


rib 




4 


9107 


1 


9/0 








9 




6 


9184 


1 


9/0 


1 


1 


8 




8 


9254 


1 


9/0 


2 


2 


7 




11 


9317 


1 


9/0 


3 


(3) 


(6) 




14 


9373 


1 


9/0 


4 


(4) 


(5) 




18 


9423 


1 


9/0 


5 


(5) 


(4) 

V / 




22 


9466 


1 


9/0 


6 


(6) 


(3) 




27/26 


9507 


1 


9/0 


7 


(7) 


(2) 




33/34 


9555 


1 


9/0 


9 








9 


9278 


2 


8/0 








8 




12 


9343 


2 


8/0 


1 


(1) 


(7) 




17 


9402 


2 


8/0 


2 








21 


9452 


2 


8/0 


3 


(3) 


(5) 




26/27 


9493 


2 


8/0 


4 


(4) 


(4) 




31/32 


9544 


2 


8/0 


5 






* 


35 


9579 


2 


8/0 


6 








41/42 


9623 


2 


8/0 


7 








19 


9428 


3 


7/0 





(0) 


(7) 




24/25 


9479 


3 


7/0 


1 


(1) 


(6) 




28/29/30 

/ / 


9513 


3 


7/0 


2 






* 


36 


9585 


3 


7/0 


3 








42/41 


9634 


3 


7/0 


4 


(4) 


(3) 




44 


9659 


3 


7/0 


5 








48/49 


9681 


3 


7/0 


6 








34/33 


9562 


4 


6/0 











43 


9655 


4 


6/0 


2 








51/50 


9717 


4 


6/0 


3 








54 


9740 


4 


6/0 


4 








81/80 


9922 


4 


6/0 


6 








46 


9665 


5 


5/0 











61 


9794 


5 


5/0 


2 








65/67 


9822 


5 


5/0 


3 








73 


9879 


5 


5/0 


4 








80/81 


9917 


5 


5/0 


5 






* 


60 


9784 


6 


4/0 











72 


9862 


6 


4/0 


1 






* 


77 


9898 


6 


4/0 


2 


(2) 


(2) 


* 


81/80 


9922 


6 


4/0 


3 






* 


86 


9967 


6 


4/0 


4 






* 


91 


10008 


6 


4/0 


5 






* 


89 


9993 


7 


3/0 


2 


(2) 


(1) 




95/96 


10047 


7 


3/0 


3 









to be continue on next page 



Number(s) 


Energy 


If = nr 
J J 


In-i-h /TIq 


t 


Ua 






97 


10053 


8 


2/0 


1 


(1) 


(1) 


* 


101 


10087 


8 


2/0 


2 








105/106 


10122 


9 


1/0 


1 


(1) 


(0) 




15 


9394 





8/1 











23 


9473 





8/1 


1 








32/31 


9547 





8/1 


2 






* 


40 


9611 





8/1 


3 






* 


47 


9670 





8/1 


4 


(4) 


(4) 




52 


9723 





8/1 


5 


(5) 


(3) 




59 


9778 





8/1 


6 


(6) 


(2) 




64 


9817 





8/1 


7 








38 


9597 


1 


7/1 











45 


9660 


1 


7/1 


1 


(1) 


(6) 




50/51 


9714 


1 


7/1 


2 








58 


9766 


1 


7/1 


3 


(3) 


(4) 




63 


9808 


1 


7/1 


4 


(4) 


(3) 


* 


69/70 


9849 


1 


7/1 


5 








74 


9883 


1 


7/1 


6 






* 


56 


9756 


2 


6/1 











62 


9802 


2 


6/1 


1 


(1) 


(5) 


* 


71 


9859 


2 


6/1 


2 






* 


78/79 


9913 


2 


6/1 


3 






* 


83 


9944 


2 


6/1 


4 






* 


87 


9972 


2 


6/1 


5 


(5) 


(1) 


* 


75/76 


9887 


3 


5/1 











90 


10002 


3 


5/1 


2 


(2) 


(3) 


* 


94 


10042 


3 


5/1 


3 


(3) 


(2) 


* 


106/105 


10125 


4 


4/1 


2 


(2) 


(2) 


* 


120 


10256 


5 


3/1 


2 


(2) 


(1) 


* 


127 


10328 


5 


3/1 


? 






* 


132/133 


10384 


6 


2/1 


? 








135/136 


10432 


7 


1/1 


? 








70/69 


9854 





6/2 









* 


92/93 


10013 





6/2 


3 








104 


10115 





6/2 


5 


(5) 


(1) 




103 


10094 


1 


5/2 


1 


(1) 


(4) 


* 


108 


10136 


1 


5/2 


2 


(2) 


(3) 


* 


113 


10177 


1 


5/2 


3 


(3) 


(1) 


* 


117 


10215 


1 


5/2 


4 






* 


118/119 


10233 


1 


5/2 


5 






* 


114 


10181 


2 


4/2 


? 






* 


147 


10588 


4 


2/2 


? 








150 


10633 


5 


1/2 


? 






* 
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Table II 



Table II: This table gives the classification and assignment of all states organized around 
fibers in b direction. The first column gives the number of the state ordered according to 
energy. If the structure only becomes clear after demixing, then it contains in addition 
the number(s) of the state(s) with which it is demixed. The second column gives the 
energy in units of cm~^. The third column gives the longitudinal quantum number Z^. 
The fourth and fifth column give the transverse quantum numbers ti and t2- 



Number(s) 


Energy 


h 


h 


t2 


161 


11096 











160 


10953 





1 





159 


10917 


1 








158/157 


10834 


1 


1 





157/158 


10818 





2 





156 


10768 


1 


2 





155 


10730 








1 


154 


10720 


1 


3 





153 


10706 





1 


1 
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Figure 1: Plots of density ( part a ) and pliase ( part b ) respectively in the plane 
V'b = V'o of density ( part c ) and phase ( part d ) respectively in the plane ijjf = 
for the state 4. In parts a and c darker grey means higher density. In parts b and d 
points with positive real part of the wave function are marked by a dot and points with 
a positive imaginary part of the wave function are marked by an open circle. All frame 
boundaries run from to 27r 
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Figure 2: Plots for state 46 in the plane tpi, = ipa + ( parts a and b ) and in plane 
ipf = Tr/2 { parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 3: Plots for state 160 in the plane ipf = —ipa ( parts a and b ) and in plane 
■i/'b = ( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 4: Plots for state 111 in the plane tpf = tph { parts a and b ) and in plane 
V'/ = 7r/2 ( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 5: Plots for state 103 in the plane ipf = { parts a and b ) and in plane tpa = 7r/2 
( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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Figure 6: Plots for state 140 in the plane ^|Ja = Sir/ 2 ( parts a and b ) and in plane 
V'b = 7r/2 ( parts c and d ). Otherwise this plot is done in analogy to Fig.l. 
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